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We have computed through order f3 21 the high temperature expansions for the nearest neighbor 
spin correlation function G(N, (5) of the classical N- vector model, with general N, on the simple-cubic 
and on the body-centered-cubic lattices. For this model, also known in quantum field theory as the 
lattice O(N) nonlinear sigma model, we have presented in previous papers extended expansions of 
the susceptibility, of its second field derivative and of the second moment of the correlation function. 

■ Here we study the internal specific energy and the specific heat C(N,/3), obtaining new estimates 
' of the critical parameters and therefore a more accurate direct test of the hyperscaling relation 

^\ , dv(N) = 2 — a(N) on a range of values of the spin dimensionality N, including TV = [the self- 

t-H ■ avoiding walk model], N = 1 [the Ising spin 1/2 model], N — 2 [the XY model], N = 3 [the classical 

' Heisenberg model]. By the newly extended series, we also compute the universal combination 

of critical amplitudes usually denoted by 7?^ (TV), in fair agreement with renormalization group 

' estimates. 

Tj- \ PACS numbers: 05.50+q, 11.15.Ha, 64.60.Cn, 75.10.Hk 

> : 

O . I. INTRODUCTION 

■ We .continue in this note the analysis of recently extended§~! high temperature (HT) expansions for the TV- vector 
modeJj with general spin dimensionality TV. Our computation is concerned with the d-dimensional bipartite lattices, 
namely the simple-cubic (sc) lattice, the body-centered-cubic (bcc) lattice and their d— dimensional generalizations. 

0^ ' In previous papers we have tabulated: i) the HT series for the zero field susceptibility x(N, P) an d for the second 
moment of the correlation function /^(TV, /3) through order f] 21 , ii) the HT series for the second field, derivative of the 
suscap,tibility %4(7V, (3) through order /3 17 , and have analysed their critical behavior in the d — 2 casea and in the^ = 3 
caseem. We have performed the computation using the (vertex-renormalized) linked cluster expansion methods and 
have produced tables of series coefficients written as explicit functions of the spin dimensionality TV with an extension 
independent of the structure and dimensionality of the lattice. More details on the derivation of the series, and on 
• • ■ the checks of validity of our results can be found in our previous papersETa. 

In this paper we examine the series expansions of the nearest neighbor correlation function G(N, 0) through order 
/3 21 , in order to update, on a range of values of the spin dimensionality TV, the direct estimates of the parameters 
describing the behavior of the specific heat C(/3, TV) on the HT side of the critical point /3 C (TV). We also update direct 
tests of the hyperscaling relation dv{N) =J2 — a{N) and estimate a relatecUmiversal combination of critical amplitudes 
introduced by Stauffer, Ferer and WortisO and later denoted by i?^(TV)0iL2l. Here a(N) is the critical exponent of 
the specific heat and v{N) is the critical exponent of the correlation length £(iV, /?). Estimates of a(N) are also 
obtained by studying the behavior of the extended series for the susceptibility x(/3, N) and for the second moment-af 
the correlation function /j,2(P,N) at the symmetrically placed anti-ferromagnetic singular point /3^(N) = — /3 c (jV)tH. 

In order to put our work into a proper perspective, it is convenient to list the HT expansions of G(N, (3) for the 
sc, the bcc and the face-centered-cuhipifcc) lattices, which were published before our extension. As well known, for 
N = 0, the N- vector model reducesOJ't^l to the self-avoiding walk (SAW) model, and the expansion of the correlation 
functioft 1 G(0, /?), simply related to the enumeration of the self-avoiding rings (or polygons) had already been computed 
in Ref.til up to order /3 19 for the sc lattice, up to order /3 15 for the bcc lattice and up to /3 13 for the fee lattice. In 
the N = 1 case, which corresponds to the spin l/2jIS|ing model, an expansion of G(l,j3) for the sc lattice has been 
obtained a few years ago by Enting and GuttmamJl3 up to order (3 21 using finit^ lattice methods. More recently, 
within the same approach, this computation has been pushed to order /3 23 in Ref£3 and then to order /3 25 in Ref.EZl. 
Also an approximate determination of the coefficient of 17 is reported in the last Reference. An expansion through 
order /3 15 for the bcc lattice, and one for the fee lattice up to order /3 12 have been tabulated in Ref.EJ. For N = 2 (the 
XY model) the available seriealj for the bcc lattice reached the order [3 11 . In the N ^-3> case (the classical Heisenberg 
model), the series for the bcc lattice, known only up to order /3 9 , is reported in RefHl 
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Finally, let us cite anjexpansion of G(N,j3), valid for general N and for all loosely packed lattices, tabulated (with 
some misprints) ip-,RcfD up to order /3 9 , which has been later extended, to models with general anisotropic pair 
interaction in RefBa. The expansion of G(N, (3) has been recently pushedB3 to order (3 15 in the case of the sc lattice, 
but no comparable effort has been devoted to the bcc lattice. In Ref.c3, an expansion to order /3 11 , valid for general 
N, had been tabulated for the fee lattice. nr~ irn 

We should finally call the readers' attention to the valuable reviews in Refs.&E3E-3 and to the accurate recalculation, 
within the Renormalization Group, (RG) approach, of the universal critical parameters of the N— vector model pet, 
formed by Guida and Zinn-JustirO. This work is based on the recently extended field theoretic expansions of Ref.Ej 
and is also accompanied by an extensive review of the available numerical and experimental data. 

The paper is organized as follows. In Sec. II we set our notation and define the quantities we shall study. 

In Sec. Ill we discuss briefly the numerical tools used for our estimates and present the results of our analysis of 
the series. These results are compared with experimental data, with earlier work on shorter HT series, with measures 
performcd-4]p-s1iOt;hastic simulations and with RG estimates, obtai ned Qj.th . pr .| b^ t he fixed dimension (FD) perturbative 
tcchniquet3'c3tll or by the Fisher- WilsorHa e-expansion approachHESE3~EilEjl2j. 

Our conclusions are briefly summarized in Sec. IV. 

The HT series expansion coefficients of the nearest neighbor correlation function G(N, (3) expressed in closed form 
as functions of the spin dimensionality N, for the sc and the bcc lattices, have been tabulated in the appendices in 
order to make each step of our work completely reproducible. For convenience of the reader, we also have explicitly 
evaluated the series coefficients for N = (the SAW model), N = 1 (the Ising spin 1/2 model), N = 2 (the XY 
model) and N = 3 (the classical Heisenberg model). 



II. DEFINITIONS AND NOTATION 



We study the TV-vector model with Hamiltonian: 

h w = -\ E (i) 

where the variable v(x) represents a iV-component classical spin of unit length at the lattice site with position vector 
x, and the sum extends to all nearest neighbor pairs of sites. 

The basic observables are the spin correlation functions. Here we shall be interested in the connected correlation 
functions (v(0) ■ v(x)) c between the spin at the origin and the spin at the site x. In particular, the nearest neighbor 
spin correlation function is defined by 

oo 

G#(N,P) = („(0) ■ v{S)) c =J24(NW- (2) 

r=0 

where 6 is a nearest neighbor lattice vector and # stands for either sc or bcc, as appropriate. 

Due to the bipartite structure of the sc and the bcc lattices, the connected correlations (v(0) ■ v(x)) c are functions 
of /3 with the same parity as the lattice distance between the spins and hence alternate expansion coefficients vanish 
identically: in particular in our expansions of G^(N,j3) to order 21 only eleven coefficients are nonvanishing. This 
is the reason why most analyses in the literature have focused on scries for the non-bipartite fee lattice which have 
no such symmetry. 

The specific internal energy is defined by 

U#(N,P) = -lG#(N,(3) (3) 

where q is the lattice coordination number. n 
If we denote the reduced inverse temperature by t#(N) = 1-/3/ flf(N), then U#(N, j3) is expected to behave asEl 

U#(N,P)^U# g (N,P)+A*(N)(r#(N)) 1 - a ^( K l + a*(N^ (4) 

when r*(N) i 0. 

As customary, in writing the asymptotic form eq. (Q) , we have explicitly indicated the presence of the non-singular 
background U&(N,f3), because the critical singularities of the specific energy are known to be generally very weak. 
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Here A^(M) denotes the critical amplitude of the specific energy, afj(N) is the amplitude of the leading singular 
correctiorEj to scaling, 0(N) is the exponent of this correction also called confluent singularity exponent. The ellipses 
represent higher order singular or analytic correction terms. Unlike the critical exponent a(N), which is universal, 
the critical amplitudes A^(N),af^(N), etc. are expected to depend on the parameters of the Hamiltonian and on 
the lattice structure, i.e. they are non-universal. Similar considerations apply to the other thermodynamic quantities 
listed below, which have different critical exponents and different critical amplitudes, but the same leading confluent 
exponent 6(N). It is known that 9(N) ~ 0.5 for small values of ML3. Having clearly indicated which quantities are 
universal, we shall often drop the generic superscript # (or its determination) in order to avoid overburdening the 
notation. Notice also that, since there is no chance of confusion, we have generally omitted the superscript + usually 
adopted in the literature for the amplitudes which characterize the high temperature side of the critical point. 

The specific heat per site, at fixed magnetic field H, is defined as the temperature derivative of the specific internal 
energy 

C H {N,P) = ±U{N,P) = 1/P±G(N,P) (5) 
where T is the temperature. As t(N) J, 0, the critical behavior of Ch(N,(3) is described by 

C H (N,j3) ~ C r H e9 (N,l3) + A c (N)(r(N))- a W (l + a c (N)(r(N)) e ^ + ....) (6) 

with A C (N) = (1 - a(N))f3 c (N)Au(N) and a c (N) = (1 + T ^4L )a u (N). Notice that our definition of a c (N) 

conforms to general usage, but differs by a factor a(N) from eq. (1.4) of RefO. 
We have also examined the susceptibility 

xWfl=5>(0).«(f)> c , (7) 

X 

the second moment of the correlation function 

/i 2 (iV,/3) =^f 2 («(0)- V (f)) c , (8) 



and the second- moment correlation length £ definedtj, in terms of x and /Z2, by 

M2(AM) 
6*(iV,/3) 



?( Ni fl = ( 9 ) 



The susceptibility x(-^ P) 1S expected to behave as 

X (N, (3) ~ A x (N)(r(N))-^ (l + a x (iV)(r(7V)) e W + ...) (10) 

as t(N) I 0. In the case of bipartite lattices x(N,(3) has also an antiferromagnetic singularity at (3f F (N) — —0 C (N), 
and, in terms of the reduced variable f(N) = 1 — (3/f3^ F (N), we should observe the energy-like behavior 

X (N,P) ~ X re g (N,P) + B x (N)(f(N)) 1 ~ a ( N ^ + ... (11) 

as f(N) I 0. 

The second moment of the correlation function is expected to behave as 

fi 2 (N,/3) ~ A^NXriN))-^- 2 ^ (l + aii (N)(r(N)) e ^ + ...) (12) 
as t(N) I 0. At the antiferromagnetic singularity, the behavior is completely similar to that of the susceptibility 

M2 (AT,/3) ~ ^(N,/3) + B,(N)(f(N)) 1 - a ^ + ... (13) 

as f(N) I 0. 

For the correlation length we have 
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C(N, /?) ~ MN)(T{N))~ vlIf) (l + a e (iV)(r(iV)) e W + ...) (14) 
as t(N) i 0, and also in this case we expect the energy-like behavior 

£(AT, /3) ~ f (JV, /?) + Bs(N)(f (N)) 1 -^ + ... (15) 



as f (AO | 0. 

The validity of the hyperscaling relation 



dv{N) = 2- a(N) (16) 



first derived by Gunton and Buckingham^ as an inequality (with the = sign replaced by >), translates into the 
universality of the amplitude combinationEl 

R+(N) = (ga(N)A c (N)y /d As(N) (17) 

where g is a geometric factor defined by g = a d /vo, with vq the volume per lattice site and a the lattice spacing. For 
the sc lattice one has g = 1, while for the bcc lattice g — 3V3/4. 

Finally, it is useful to recallH that, in the large N limit, R^(N) » (j^) 1 ^ 3 and that Bervillier and GodrechefEl 
proposed a simple approximate extension of this relationship to small nonzero values of N in the form R^(N) w 

"W(&) 1/B - 



III. COMMENTS ON THE ANALYSIS OF THE SERIES 



A. Estimates of the specific heat exponents 



The main difficulty in computing the specific heat exponents is that a(N) is small for AT < 1 and it becomes 
negative for N > 2. Therefore the specific heat is very weakly divergent for N < I, whereas it has only a finite cusp 
for N > 2. The simplest Pade approximant (PA) techniques for estimating the critical parameters are thus expected 
to be inefficient in the former cases and completely inadequate in the latter. Moreover, it is not particularly helpful to 
differentiate the present specific heat series with respect to /3 in order to sharpen the singularity, because the extrap- 
olations become more sensitive to non-asymptotic or confluent singularity effects. In principle, the inhomogeneous 
differential approximants (DA) (thoroughly described in Refs.c2l) should perform much better than the PA's since they 
are able to detect even weak singularities and might allow, to some extent, for the confluent corrections to scaling. 
However, even after our extension of the HT series, the nonzero expansion coefficients are not sufficiently many that 
these numerical tools can be used effectively. In order to improve the precision of our estimates, we have mainly used 
simple first order DA's and have biased them—with the critical temper atvLres reliably known from our previous study 
of the strongly divergent susceptibility serieso or from other sourcesEJc3. In the particular case of the sc spin ]J2 
Ising model, we have taken advantage in our analysis also of the two additional series coefficients provided by Ref.ti 

An accurate measure of the scaling correction amplitudes of the specific heat presently seems beyond reach, although 
their qualitative behavior as functions of N is clear and completely analogous to that of a x (N) and of a^(N). More 
precisely, a^(AT) is small and negatiiasJjor N < 2, while it is positive and increasing for N > 2. Let us recall that, for 
small values of AT, RG computationsoEi indicate that the universal ratios ac(N)/a x (N) and ac ( AT) /a^ ( N) are of the 
order of the unity. On the other hand, our HT analysis of x(/3, N) and £(/3, N) suggested that a x (N) and a^(N) are 
small (negative for A^ < 2 and positive otherwise), therefore it is reasonable to neglect the corrections to scaling at the 
present level of accuracy in the specific heat series analysis. We also recall that it was convincingly inferredpLn-Ref.H 
that ac is negative in the sc, bcc and fee spin 1/2 Ising models and, in the sc case, it was suggested in Refs.E3il3 that 
ac is very small. 

Our direct estimates of a(N) from the specific heat series for the sc and the bcc lattices have been reported in 
Table 1. We have also included in this Table the values of a(N) obtained by studying the energy-like behavior of 
the susceptibility eq.(pd|) at the antiferromagnetic singularity. The study of the second correlation moment eq.jl^) 
does not produce results of comparable quality. In this computation, we have found most convenient to analyse the 
derivative of \ by second order DA's biased with the singularities at C (N) and f3^(N). Although the expansion of \ 
is effectively longer than that of the specific heat, it is not easier to measure accurately the exponent of the very weak 
antiferromagnetic singularity. Therefore the estimates of a(N) so obtained are consistent with, but not more accurate 
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than the others. In particular, we agree with the earlier estimates a(l) y ,0,105(7) and a(l) = 0.11(2) obtained by 
studies of the susceptibility for the teifift model on the bcc lattice in Refs.liilc3. 

In recent studies of the N = 1 caset^O, it has been suggested that the behavior of the specific heat series coefficients 
as functions of their order is sufficiently smooth that the traditional (biased) ratio techniques can be practically as 
accurate as the DA procedures. This remains true only for not too large values of N, since an asymptotic regime 
seems to set in later for larger N. Moreover, for N > 4, the ratio sequences show an increasing curvature indicating 
that the confluent corrections to scaling cannot be neglected anymore and therefore longer series are needed for a 
reliable analysis. I— ■ 

We have used the simplest ratio formulas, since the more elaborate variants proposed in Ref.Ea do not presently 
make much difference. If we set Ch(N, (3) = Y^=i c n(N)f3 2n , and allow for the dominant corrections to scaling with 
exponent 9(N), the ratio of the successive coefficients of the specific heat expansion in powers of (3 2 is expected to 
behave as 

= /3 2(i + I^ + Jl_ + 0( J_)) ( i8) 



Cn+1 

Therefore a can be estimated from the sequence 



= l-(^-l)n = a+^ + 0(-) (19) 



The extrapolation of these estimators to n — > oo is the main difficulty with this procedure. For N < 4, the 
estimators, when plotted versus 1/n, show only a small curvature. Therefore we have neglected the scaling correction 
b/n e and have simply taken the linear extrapolant na n — (n — l)a„_i of the last two estimators as our final estimate 
of a(N). We have then assigned very conservative uncertainties to these results (also allowing for the errors in (3 C (N)) 
and, for N > 2, we have indicated by asymmetric errors the effects of some curvature in the estimator plots. 

In Table I, we have also included the results of a few recent direct studies of the specific heat by stochastic methods. 
These studies are subject to difficulties analogous to those met in HT analyses. As a consequence, for instance, the 
MonteCarlo (MC) determinations of a(0) on the sc lattice is approximately three standard deviations away from 
the other quoted values. (We have quoted the sum of the systematic and the statistical errors separately reported in 

Rcf.y.) 

Also the value of a(l) emerging from a most accurate (see Ref.El) MC study of the sc lattice Ising model performed 
by a dedicated processor, shows a considerable uncertainty. The central value, but not the error, is somewhat improved 
(a(l) = 0.113 ± 0.023) by turning to a particular spin 1 Ising model designed to have small corrections to scaling. 

For N > 2, it is even harder to determine the exponent a(N) in MC simulations, because of the ambiguity in-thc 
separation of the non-divergent singular part of the specific heat from the regular background, ras argued in Refo. 

We have also reported a few experimental measurements of the specific heat exponentll3~E3 available for TV = 
1,2,3,4. 

In order to show quantitatively the validity of the hyperscaling relation eq. (|l6|) , our direct estimates of a(N) have 
been compared with the quantity 2-,— 3v(N) also reported in Table I and computed either from our extended HT 
expansions of the correlation lengtto for the sc and the bcc lattices, or from the estimates-Xiili^TV) obtained in the 
RG approach by fifth order e-expansion and by seventh order FD perturbation expansiono'EEI. In conclusion, the 
hyperscaling relation dv{N) = 2 — a(N) appears to be reasonably well verified within the uncertainties of the data. 



B. Estimates of R+ (N) 



We have computed the hyperuniversal combination of critical amplitudes R^(N) by two methods. In the first 
procedure, we evaluate the HT expansion of the quantity 



F 



W ffl-w(«^'(sn£)-^- « 

R+(N) 3 t(N) 2 ^ n ^ 3 ^ n \i + 0(t(N) 6 ^)) 

at the critical temperature. This computation also provides a good test of hyperscaling: indeed F(N, (3 C (N)) = 
Rf(N) 3 , if eq.© holds. Here we have found convenient to use the "simplified" first order DA's, biased with (3f(N) 
and 0(N), as described in Ref.i, and have taken the estimates of 8(N) from RefR We have reported in Table II 
only the results obtained by this method which is very stable and seems to be fairly accurate. In this case, our error 
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estimates have to allow only for the spread of the approximants as well as for the uncertainties of (3 C (N), v(N) and 
0(N). The errors quoted mainly derive from the uncertainties in 0(N), assumed to be generally of the order of 10% 
and from the uncertainties of v(N). The estimates of R^(N) obtained by PA's of F(N, /3 C (N)) are systematically 
smaller by « 5%, indicating, in our opinion, that the "simplified" DA's are likely to allow more accurately for the 
sizable negative amplitude corrections to scaling. The usual first order DA's biased with f3 c (N) also seem to lead to 
less accurate estimates. 

In the second approach, we obtain R^(N) from eq.(17), after computing separately Ac(N) and A^(N) from 
the specific heat and the correlation length series respectively, by DA's biased with the critical temperatures and 
exponents. This second method leads to results systematically smaller (by »1- 2%), than those reported in Table 
II and it is subject to a larger uncertainty, due to the necessity of biasing the direct computation of Ac(N) also with 
the exponents a(N), whose relative error may be considerable. n 

In the same Table we have also reported the values of R^(Nl computed via RGu either to second order in the 
e-expanaiantLor to fifth order in the FD perturbation expansionc3. We have also included earlier estimates obtained 
in Refs.Biij'EIJ from the analysis of shorter HT series, by the second above mentioned method. 

A recent MC simulation!! of the Ising model on the sc lattice has determined the universal quantity /f c (l) ^A| c (l) 
which is closely related to R^(l). Here /| c (l) denotes the amplitude of the singular part of the free energy. For 
convenience, we have translated this result into the estimate of -R^(l) reported in Table II, by using the value 
a(l) = 0.1076(30), obtained in the same Ref.S from the hyperscaling eq. (|l^) . 

The values from the approximate formula of Bervillier and Godreche have been obtained assuming for v(N) the 
FD perturbative results of Ref.Ej. We are unable to give sensible error estimates in this case, but it interesting to 
quote at least the uncertainties deriving from those of v(N). 

Finally, we should mention that, to our knowledge, no other evaluations of R^(N) for N — and N = 4 are quoted 
in the literature. 



C. Estimates of non- universal critical parameters 

In Table III, we have reported our estimates of some non-universal critical parameters, for various values of N. 
The inverse critical temperatures Pf(N\n which have been always used in the biased analyses of this paper were 
determined in Reffl or taken from RefsJidH. n 

The critical amplitudes A^(N) of the second-moment correlation length were determined in Ref.B. 

The critical specific energies U&(N, (3 C ) and the critical values of the regular part of the specific heat Cf eg (N, f5 c {N)) 
have beeta-, obtained by first order DA's biased with (3 C (N). Also these data are compatible with the earlier 
estimatesE3£il. 

We have computed the critical amplitudes of the specific heat Af^(N) in two ways: either indirectly, namely from 
our estimates of R^(N) by using the knowledge of Ajt(N) and of a(N), or directly, from the specific heat by DA's 
biased with f3 c (N) and a{N). The two methods yield compatible results. We have chosen to report in Table III the 
results of the first approach. Therefore the relatively large errors of Aq(N) mainly reflect the uncertainty of a(N), 
which, for N = 2, is so considerable that it is not useful to report any estimates in this case. (For the same reason we 
have not reported estimates of Cf eg {2, (3 C {2)).) On the other hand, the uncertainties of the products a(N)A*{N) are 
more modest and therefore it can be of some interest to quote our estimates for N = 2, namely a(2)A^f(2) — 0.42(1) 
and a(2)A<g c (2) = 0.44(1). 

We should stress that here the meaning of the errors for Rf(N), A%(N) etc. is not the same as in earlier studies, 
where the errors describe the spread of the estimates in computations performed at sharply fixed values of a(N) and 
ft c (N). If, in those computations, we allowed also for the uncertainty of a(N), then the estimates and the errors 
of R^(N), Afj(N) etc. would become completely compatible with our results. Therefore, for instance, we have 

reported in Table III the central values of the estimates of A*(l) from Ref.0, based on the assignments a(l) = 0.104, 
/3* c (l) = 0.221630 and I3 h c cc {l) = 0.157368, but we have taken the liberty of suggesting much lar ger errors, which 
correspond to an indicative 5% uncertainty of a(l). 

Finally, it is interesting to notice that the product a(N)A^(N), which is derived with good accuracy from R^(N), 
remains positive in the range of N examined here. Therefore, when a(N) changes its sign for N = N < 2, the same 
must happen for Aq(N). Analogously Cf eg {N, /3 C ), which is negative for N = 0, 1, has to change sign for N > N, in 
order that the maximum of the specific heat stays positive. 
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IV. CONCLUSIONS 



We have analysed our extended HT expansion of G(N, f3) for the sc and bcc lattices in order to update the direct 
estimates of the critical exponent a(N) and of the hyperuniversal combination of amplitudes Rf(N), over a range of 
values of N. 

Due to the smallness of a(N) and to the limited effective length of the series, the relative accuracy of our ex- 
trapolations is still generally inferior to that already achieved in our recent HT studies of the susceptibility and of 
the correlation length critical exponents^. However, within the error limits, the main predictions of universality, 
hyperuniversality and hyperscaling appear to be well verified and the overall agreement between the HT and the RG 
estimates of the universal observables is good. 
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TABLE I. In the first six lines we have reported the direct HT estimates of the critical exponents a(N) obtained in this 
work by various routes: by first order DA's of the specific heat biased with f) c (N); by similarly biased extrapolation of ratios of 
the specific heat series coefficients and by second order DA's of dx/df3 biased with p c (N) and /3" (N). We have then reported 
earlier direct estimates from shorter HT series, some direct MC determinations, and a few experimental measures. For each 
value of N, our estimates of a(N) have to be compared with the quantity 2-3v(N) reported in the last four lines and obtained 
either from our previous HT study of the correlation length series or from RG estimates via e-expansion and via FD perturbative 
expansion. 



N 







1 



Cg{N,0) DA 
C^ c (N,/3) DA 
C|f C (Ar,/3)Ratio Ext. 
C!T(iV,/3)Ratio Ext. 
%-{N,P) DA 
*$P-(N,P) DA 
X 6cc (l,/3) DA Ref. 

C H {N,/3 

MC Refs^ YlHl 
Exper. RefsBO 
2-3v(N) DA sc Reft 
2-3v(N) DA bcc Ref.0 , 
2-3f (iV) e-expans. Ref* 
2-3u(N) FD-pert. Ref.E 



,cc (l,/3) DA Ref. 

f. 



0.24(1) 

0.23(1) 

0.236(8) 

0.233(8) 

0.239(8) 

0.233(6) 



0.275(15) 

0.2366(18) 
0.2363(18) 
0.2375(54) 
0.235(3) 



0.103(8) 

0.105(9) 

0.104(6) 

0.106(6) 

0.13(3) 

0.107(8) 

0.105(7) 

0.105(7) 

0.125(25) 

0.125(23) 

0.107-0.112 

0.1055(24) 

0.1076(15) 

0.1121(78) 

0.109(4) 



-0.014(9) 

-0.019(8) 

-0.020(8) 

-0.022(6) 

0.02(3) 

-0.01(2) 



-0.02(3) 

-0.01285(38) 

-0.025(6) 

-0.022(6) 

-0.0055(120) 

-0.011(4) 



-0.11(2) 

-0.13(2) 

-0.15(±J) 

-0.16(+?) 

-0.13(3) 

-0.138(8) 



-0.22(4) 

-0.23(16) 

-0.135(2) 

-0.148(6) 

-0.142(6) 

-0.115(16) 

-0.122(10) 



+t l) 



-0.22(4) 
-0.25(3) 
-0.27( 
-0.29( 
-0.24(3) 
-0.23(2) 



-0.20(5) 

-0.277(9) 

-0.268(9) 

-0.211(24) 

-0.223(18) 



TABLE II. Estimates of the hyperuniversal quantity R£(N). The results of our HT series computation are compared with 
RG estimates via e-expansion or via fixed-dimension perturbative expansion, with a heuristic approximate formula and with 
experimental measures. 



N 

HT sc (this work) 
HT bcc (this work), 
RG e-expans. Ref.B,|— ■ 
RG FD-ij^URefsBEl 
HT RefeMBB 



MC Ref. 
Exper. Refs 




_0 

0.258(3) 
0.258(3) 



0.273(4) 


0.361(4) 


0.431(5) 


0.272(4) 


0.362(4) 


0.433(5) 


0.27 


0.36 


0.42 


0.270(1) 


0.361(2) 


0.435(2) 


0.2659(7) 


0.36(1) 


0.42 


0.2685(60) 






0.271(1) 


0.363(1) 


0.439(2) 


0.25-0.32 




0.40-0.45 



0.497(6) 
0.500(6) 



0.506(4) 
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TABLE III. Estimates of non-universal parameters. We report the critical inverse temperatures f3f(N) always used in our 
biased procedures, the critical amplitudes A^(N) and A^(N), the critical specific energies U*(N,f3 c ) and the critical values 
of the regular part of C^(N, (3). 

N 1 2 3 4 

WW) HT Rcfs %*i 

/3 bcc {N) HT RefrijE] 



Ai c (N) HT Ref.t 
l ? "(N) HT Ref. 



A°S{N) (this work) 



1 C 



: (JV) (th 



A%(N) MC Ref « 
A%(N) HT Ref.t 
A b S c {N) HT Ref. f 11 
U 3C (N,p c 
U bcc (N,f3, 
U sc (N,f3 c 
U ac (N,f3 c 



(this work) 
(this wipk) 
HT RefM 
MC RefM 
HT Ref.B 
i sc (this work 
i MC sc Refs, 
I MC sc Ref.t 
C™ 9 {N,f3 c ) bcc (this work) 
C r * g (N,p c ) MC bcc RefM 



C r H ea (N,(3, 
C r H e9 (N,f3, 



(N,0 C 



213493C3") 


2216544(31 


45419(31 


69305(41 


0.153128(3) 


0.157373(2) 


0.320427(3) 


0.486820(4) 


5101(31 


5027(31 


4814(31 


4541(31 


4846(21 


4659(21 


4371(21 


4072(21 


546(81 


1.49(5) 




-6.0(6) 


481 (61 


1 43(41 




-6.5(6) 




1.45(9) 








1 464(901 








1 431 (801 






-1.004(3) 


-0.991(1) 


-0.990(3) 


-0.991(3) 


-1.0990(2) 


-1.0903(6) 


-1.0896(8) 


-1.0919(4) 




-0.9902(1) 








-0.9904(8) 








-1.0904(1) 






-0.66(3) 


-1.67(3) 




4.9(4) 




-1.64(11) 




5.79(12) 








5.70(12) 


-0.68(2) 


-1.64(3) 




5.2(4) 



0.93600(4) 
0.65542(3) 
0.4155(3) 
0.3691(2) 
-6.5(3) 
-7.2(3) 



-0.994(4) 
-1.0951(2) 



4.2(3) 
4.3(3) 



5.54(14) 
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APPENDIX A: THE NEAREST NEIGHBOR CORRELATION FUNCTION ON THE SC LATTICE 

The HT expansion coefficients of the nearest neighbor correlation function on the sc lattice are 

K ' TV 
8 + 3TV 

a 3 (N) 



N 3 (2 + TV) 



352 + 168JV + 22TV 2 

U5<N) = 

TV 5 (2 + TV)(4 + TV) 

_ 105984 + 1547527V + 85056TV 2 + 21960TV 3 + 2754TV 4 + 135TV 5 
a7( ' ~ TV 7 (2 + TV) 3 (4 + 7V)(6 + TV) 

For the coefficients which follow, it is typographically more convenient to set a r (N) — P r (N)/Q r (N) and to tabulate 
separately the numerator polynomial P r (N) and the denominator polynomial Q r (N), 

P 9 (7V) = 12349440 + 178713607V + 100102407V 2 + 27516807V 3 + 4057767V 4 + 308767V 6 + 9547V 6 
Q 9 (7V) = 7V 9 (2 + 7V) 3 (4 + 7V)(6 + 7V)(8 + TV) 
P lx (TV) = 124861808640 + 364560318464TV + 467027804160TV 2 + 345395589120TV 3 + 163465120768TV 4 + 5 19376629767V 5 + 



Qn(7V) = N X1 (2 + TV)°(4 + TV) 3 (6 + TV)(8 + TV)(10 + TV) 
P 13 (N) = 24917940633600 + 71794651299840TV + 91099400634368TV 2 + 67066306363392TV 3 + 31821500096512TV 4 + 
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102421285908487V" + 22953204715527V" + 3617895637767V' + 399248565127V° + 30149464647V + 
1480813127V 10 + 42497127V 11 + 53892JV 12 
Ois(JV) = 7V 13 (2 + 7V) 5 (4 + 7V) 3 (6 + 7V)(8 + TV) (10 + N)(12 + TV) 
(TV) - 867654721119191040 + 3616829986427633664TV + 68915837394286018567V 2 + 79573832548378214407V 3 + 
6225913571872604160Af 4 + 3498334649912000512 N 5 + 1460523056889888768 W 6 + 
4625632235925667847V 7 + 1125211548203499527V 8 + 211542535316848647V 9 + 
30762403605872647V 10 + 3443764911744007V 11 + 293392594145607V 12 + 18634096654567V 13 + 
852237782567V 14 + 26444517687V 15 + 496791147V 16 + 4250077V 17 
Qib(JV) = 7V 15 (2 + 7V) 7 (4 + 7V) 3 (6 + 7V) 3 (8 + 7V)(10 + 7V)(12 + 7V)(14 + TV) 
P 17 (7V) = 3948322260048528015360 + 182265982596873254338567V + 389880217237899360337927V 2 + 
513238696901276451471367V 3 + 465835507424588338298887V 4 + 309606814623706518650887V 5 + 
156232516353352794112007V 6 + 61261147711923599441927V 7 + 18951343400756279377927V 8 + 
4670228089812312227847V 9 + 921861818644423516167V 10 + 146036835964908257287V 11 + 
18538630987151370247V 12 + 1876060642026608647V 13 + 149886695254955527V 14 + 
9308100122144647V 15 + 438623233288647V 16 + 15105378825927V 17 + 
357260754727V 18 + 5165868767V 19 + 34266107V 20 
Qi 7 (7V) = 7V 17 (2 + 7V) 7 (4 + 7V) 5 (6 + 7V) 3 (8 + 7V)(10 + 7V)(12 + 7V)(14 + 7V)(16 + TV) 
P 19 (7V) = 330768394077031316324352000 + 19214894928068294618382336007V+ 
52443527485608930541200998407V 2 + 89435854981410476078923776007V 3 + 
106928579326845764041388851207V 4 + 95337461125086677039222620167V 5 + 
65840534257305886001998069767V 6 + 36116025803779273901735936007V 7 + 
16012077391462158006988308487V 8 + 5808922314056280184308367367V 9 + 
1739582312375680987491205127V 10 + 432634818212640258592604167V 11 + 
89701597975609364619591687V 12 + 15535633987901684283146247V 13 + 
2247894974205115799633927V 14 + 271277170914385267343367V 15 + 
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2720424198488158732288JV 16 + 2253266792765734184967V 17 + 152765305959025853447V 18 + 
8368844711097224967V 19 + 363715847042973447V 20 + 12212318316035527V 21 + 
304415355645767V 22 + 5282113357527V 23 + 56660577527V 24 + 281133667V 25 
Q 19 (7V) = TV 19 (2 + 7V) 9 (4 + 7V) 5 (6 + 7V) 3 (8 + N) 3 (W + 7V)(12 + 7V)(14 + 7V)(16 + iV)(18 + AT) 
P 2i (7V) = 122469323387965953278371430400 + 7039681357139968749683186073607V + 
19020187744014823729253811486727V 2 + 32128331430094756215860011991047V' 3 + 
38074215457815581892635661434887V 4 + 33677732917117004125208313856007V 5 + 
23098565275274440830011404451847V 6 + 12600258116295110348045560053767V 7 + 
5564347967448682315291515944967V 8 + 2014626193173915173739099586567V 9 + 
603523336261384848497537187847V 10 + 150572129813681896217892618247V 11 + 
31424582320403625319623557127V 12 + 5500720004891682442809507847V 13 + 
808419835958253762373058567V 14 + 99692537424225810364743687V 15 + 
10291774527135660977479687V 16 + 885654567153173906063367V 17 + 
63111616662732417105927V 18 + 3688584807784925132807V 19 + 
174435654196059112967V 20 + 6548519215100171527V 21 + 189874468392175367V 22 
+4080991742348487V 23 + 60858408716807V 24 + 557730637927V 25 + 2339665567V 26 
Q 21 (7V) = N 21 (2 + 7V) 9 (4 + 7V) 5 (6 + 7V)' 3 (8 + 7V) 3 (10 + 7V)(12 + 7V)(14 + 7V)(16 + 7V)(18 + N)(20 + TV) 

In particular, for N = [the SAW model], we have (in terms of the variable = (i/N) 

G(0, 0) = + 4/3 3 + 44/3 5 + 552/3 7 + 8040/3 9 + 127016/3 11 + 2112320/3 13 + 36484128/3 15 + 648529392/3 17 
+ 11790401800/5 19 + 218273957968/3 21 + ... 

For N = 1 [ the spin 1/2 Ising model], we have 

G(l, 0) = + 11/3/3 3 + 542/15/3 5 + 123 5 47/ 3 1 5/3 7 + 144 734 42 / 2 8 35/3 9 + 11336607022/155925/3 11 + 605337636044/552825,3 13 + 
10976336338579019/638512875/3 15 + 3022947654230404442/10854718875/3 17 + 8582760723898537620322/1856156927625/3 19 + 

15262009695163033631128084/194896477400625/3 21 + ... 

For N = 2 [ the XY model] , we have 

G(2,0) = 1/2/3 + 7/16/3 3 + 97/96,3 5 + 5103/2048/3 7 + 459719/61440/3 9 + 218788559/8847360/3 11 + 3579816967/41287680,3 13 + 
20154248931151/63417876480/3 15 + 4126827327908711/3424565329920/3 17 + 2142771095208749011/456608710656000,3 19 + 

562665453010146198199/30136174903296000/3 21 + ... 

For N = 3 [the classical Heisenberg model], we have 

G(3,0) = 1/30 + 17/135/3 3 + 1054/8505/3 5 + 80909/637875/3 7 + 95738/601425/3 9 + 5992817408726/27152760009375£J 11 + 

11357358327572/34910691440625/3 13 + 156550175755271443/311577921107578125/3 15 + 
190956190202826883834/237337400087308828125/3 17 + 56535690823720347706912558/42645970734688086781640625/3 10 + 
358752594209204675460504716/160503926219644253887265625/3 21 + ... 
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APPENDIX B: THE NEAREST NEIGHBOR CORRELATION FUNCTION ON THE BCC LATTICE 

The HT expansion coefficients of the nearest neighbor correlation function on the bcc lattice are 

24 + 117V 

a 3 (N) 



a 5 (JV) 



7V 3 (2 + TV) 

1776 + 1044TV + 1527V 2 
TV 5 (2 + TV)(4 + TV) 



_ 1050624 + 1713024TV + 1062432TV 2 + 312600TV 3 + 44090TV 4 + 2395TV 5 
a7(Ar) ~ TV 7 (2 + TV) 3 (4 + N)(6 + TV) 

For the coefficients which follow, it is typographically more convenient to set a r (N) = P r (N)/Q r (N) and to tabulate 
separately the numerator polynomial P r (N) and the denominator polynomial Q r (N), 

P 9 (TV) = 237680640 + 391630080TV + 251136960IV 2 + 79995360TV 3 + 13572456TV 4 + 1175956TV 5 + 40904AT 6 
Q a (N) = N 9 (2 + Nf(4 + JV)(6 + N)(8 + N) 
Pu(N) = 4657615994880 + 146624394362887V + 20306810757120Af 2 + 16297064577024Af 3 + 8408736450048JV 4 + 
2927305709568W 5 + 701958299776Af 6 + 116098602304JV 7 + 130014820807V 8 + 9405463047V 9 + 396188967V 10 + 7371127V 11 

Qn(7V) = AT 1X (2 + 7V) 5 (4 + 7V) 3 (6 + 7V)(8 + 7V)(10 + TV) 
P 13 (7V) = 1804392176025600 + 56604209047142407V + 78383698931220487V 2 + 63201160293089287V 3 + 32991742874173447V 4 + 
11738998724060167V 5 + 2921019880949767V 6 + 512946695786887V 7 + 63227376982727V 8 + 
5347494982887V 9 + 295206408087V 10 + 9569574407V 11 + 137992327V 12 
Ois(JV) = 7V 13 (2 + 7V) 5 (4 + 7V) 3 (6 + 7V)(8 + 7V)(10 + 7V)(12 + TV) 
P 15 (7V) = 122002510248369192960 + 543062014542747795456TV + 1106788272284626845696TV 2 + 
1369759313298192334848TV 3 + 1151523649799700086784TV 4 + 697153536634263011328TV 5 + 



58543844261560627207V 9 + 930798833517987840TV 10 + 114224602657806848TV 11 + 
106968880312488007V 12 + 7491476163933287V 13 + 379272136091687V 14 + 
13091428536247V 15 + 275304441147V 16 + 2657766997V 17 
Q 15 (7V) = 7V 15 (2 + 7V) 7 (4 + 7V) 3 (6 + 7V) 3 (8 + 7V)(10 + 7V)(12 + 7V)(14 + TV) 
P 17 (7V) = 1078176657766748635791360 + 5311087820388289065517056TV+ 
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12134778449671899486093312AT 2 + 17084121119131816471560192iV 3 + 
16608988472236389552881664AT 4 + 11844719796924403830226944JV 5 + 
6426005983019854016544768AT 6 + 2714739161543420662382592AT 7 + 
906815315440850230886400A rS + 241845902161952176398336JV fl + 
51782172208072962975744JV 10 + 8918298897939850518528^" + 
1233687735206166823424JV 12 + 136374954925264681472JV 13 + 
11933805821587623936JV 14 + 814292868353822272A rls + 42326351689562720JV 16 + 
1615854325367776JV 17 + 42636630380712JV 18 + 693565371332JV 19 + 5232689960JV 20 
Q 17 (N) = N 17 (2 + N) 7 (i + Nf(6 + N) 3 (8 + N)(W + N)(12 + AT) (14 + AT)(16 + N) 
P 19 (N) = 175433838338452762972599091200 + 1078203808265217149807540305920IV+ 
3116173607652516436321480212480JV 2 + 5633274990887456528054857236480A' 3 + 
7147764845163022335548483174400JV 4 + 6772123960394212475971911548928A r5 + 
4976697709965417132847811002368JV 6 + 2909196906278057822968464015360Af 7 + 
1376618167955575605597824876544JV 8 + 533877242182341657037822230528 W 9 + 
171186415133933732406244147200A' 10 + 45658929843087077741062520832JV 11 + 
10169112743416038967608705024AT 12 + 1894893729303875646595006464A r13 + 
295461155711988902515834880JV 14 + 38487351224017884772339712Af ls + 
4173283404716049615550464JV 16 + 374483776338165994216448A rl7 + 
27568453015797659725824JV 18 + 1644491785569570191872Af 19 + 
78098433373280888576JV 20 + 2878619676002249280JV 21 + 
79247785021379008iV 22 + 1531012629840624JV 23 + 18487394459632AT 24 + 104841714952AT 25 
Q 19 (N) = JV 19 (2 + Af) 9 (4 + JV) 5 (6 + AT) 3 (8 + JV) 3 (10 + N)(12 + AT) (14 + N)(W + N)(IS + N) 
P 21 (N) = (126175500888039348819208018329600 + 771182011973845295293568951255040^+ 
2217488378116952516842355989413888AT 2 + 3990433665646539615566917678399488AT 3 + 
5043598010921343492445826797535232AT 4 + 4763919940749294277936547050291200AT 5 + 
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3493715661682982610876998803783680A r<i + 2040596945264227996384033763229696 N 7 + 
966223362153625194801269330411520Af 8 + 375629302981996871057940890517504AT 9 + 
120996637012851297773150580768768Af 10 + 325040728550728192742675707330567V 11 + 
7314027149319445947024024403968A f12 + 1382157562600983080851995279360A rl3 + 
219565426180384061435124916224Af 14 + 29302944828757838601074274304AT 15 + 
3278043851101943224116459520Af 16 + 306114536101034343439527936iV 17 + 
23712149887809097730532352AT 18 + 1509817181912122244980224JV 19 + 
78017870011349442092288AT 20 + 3213483265751133834688Ar 21 + 102816420657321623712Af 22 + 
2458334158401005552Af 23 + 41257373964220632AT 24 + 432739125346952AT 25 + 2130772922816A f26 
Q 2 i(N) = N 21 (2 + N) a (4 + Nf(6 + N) 3 (8 + N) 3 (W + JV)(12 + JV)(14 + JV)(16 + JV)(18 + W)(20 + N) 

In particular, for N = [the SAW model], we have (in terms of the variable (3 = (3/N) 

G(0, 0) = + 12/3 3 + 222/3 5 + 5472/3 7 + 154740/3 9 + 4737972/3 11 + 152960220,3 13 + 5130099672,3 15 + 177095284092^" 

+6253425298080/3 19 + 224879383796232/3 21 + ... 

For N = 1 [ the spin 1/2 Ising model], we have 

G(l, 0) = + 35/3/3 3 + 2972/15/3 5 + 279011/63/3 7 + 46439636/405/3 9 + 100877055128/31185/3 11 + 587703506650264/6081075/3 13 + 
10981652882712713/3648645/3 15 + 1049923978894758374012/10854718875/3 17 + 1182698210781462071363672/371231385525/3 19 + 

2980059927747623321534851312/27842353914375/3 21 + ... 

For N = 2 [ the XY model], we have 

G{2,0) = 1/2/3 + 23/16/3 3 + 559/96/3 s + 187645/6144/3 7 + 11417419/61440/3 9 + 10934199853/8847360/3 11 + 
1081218105839/123863040/3 13 + 4085878131871327/63417876480/3 15 + 1683908448367350071/3424565329920/3 17 + 
753925204192677068291/195689447424000/3 19 + 929152049503798552678997/30136174903296000/3 21 + ... 

For N — 3 [ the classical Heisenberg model], we have 

G(3,/3) = 1/3/3 + 19/45/3 3 + 2092/2835/3 6 + 349939/212625/3 7 + 2147444/505197/3 9 + 108732988464808/9050920003125/3 11 + 

9339742669288/258597714375/3 13 + 35412600932786885263/311577921107578125/3 15 + 
35816645375345371477924/96693014850385078125/3 17 + 5880568448944900843943527784/4738441192743120753515625/3 19 + 
1137371495914136837811604445344/267506543699407089812109375/3 21 + ... 
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